Nonlocal effects in high energy charged particle beams 
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Within the framework of the thermal wave model, an investigation is made of the longitudinal 
dynamics of high energy charged particle beams. The model includes the self-consistent interaction 
between the beam and its surroundings in terms of a nonlinear coupling impedance, and when 
resistive as well as reactive parts are included, the evolution equation becomes a generalised nonlinear 
Schrodinger equation including a nonlocal nonlinear term. The consequences of the resistive part 
on the propagation of particle bunches are examined using analytical as well as numerical methods. 
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I. INTRODUCTION 

The thermal wave model (TWM), 0,110], describes 
the dynamics of high energy charged particle beams in 
accelerators. In the TWM approach, the beam is charac- 
terised by a complex-valued wave function, which sat- 
isfies a Schrodinger-like evolution equation, where the 
beam emittance replaces Planck's constant, and the in- 
tensity of the wave function corresponds to the beam par- 
ticle density. The Schrodinger potential, which describes 
the interaction between the beam and its surroundings, 
can be expressed in terms of a coupling impedance, and 
due to collective effects, the coupling is a nonlinear func- 
tion of the beam density. For purely reactive impedances, 
the TWM equation reduces to the well-known nonlinear 
Schrodinger equation. However, by including also the re- 
sistive part, the evolution equation becomes a generalised 
Schrodinger equation containing a new term, which is 
both nonlinear and nonlocal. The modulational insta- 
bility properties of this new equation have been analysed 
previously, j^j , and have been shown to agree with results 
obtained using classical approaches, including kinetic ef- 
fects like Landau damping, 5]. 

In the present work we consider the longitudinal dy- 
namics of particle bunches under the influence of the cou- 
pling impedance. Since the case with a purely reactive 
impedance is well known, main emphasis is on the situa- 
tion where the resistive part is included. The dynamical 
evolution then proceeds as a competition between linear 
diffraction, nonlinear self-focusing/defocusing, and non- 
local self-steepening. It is found that the bunch is accel- 
erated/decelerated, and the self-steepening effect makes 
the pulse shape asymmetric with an extended tail, and 



eventually a wave-breaking phenomenon can appear on 
the steepening edge. An approximate solution for the 
dynamics of soliton-shaped bunches in the presence of 
a small resistive impedance is found using variational 
methods, and a perturbation solution for the initial evo- 
lution of a particle bunch is given. This illustrates the 
interplay between the different effects and is qualitatively 
in agreement with results obtained by other means, as 
summarised, e.g., in 0. No stationary solutions with fi- 
nite particle number are possible, but semi-infinite shocks 
and pulses with extended wake fields are found. The an- 
alytical predictions are confirmed by numerical simula- 
tions of the full generalised nonlinear Schrodinger equa- 
tion. 



II. THE GENERALISED NONLINEAR 
SCHRODINGER EQUATION 

Within the TWM, the longitudinal dynamics of par- 
ticle bunches are analysed in terms of a complex beam 
wave function tp{x,z), where z is the distance of propa- 
gation and X is the longitudinal extension of the particle 
beam, measured in the moving frame of reference. The 
particle density, n{x,z), is related to the wave function 
according to n{x,z) — \^{x,z)\^, see 0. The evolution 
of the beam is described by the generalised Schrodinger 
equation 
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where the longitudinal diffraction parameter, a, can be 
both positive and negative, depending on the phase slip 
parameter, and the nonlinear parameters, k and /i, 
are proportional to the imaginary (reactive) and real (re- 
sistive) parts, respectively, of the coupling impedance. 
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In the case ^ = 0, Eq. reduces to the fundamen- 
tal nonhnear Schrodinger equation for which a wealth 
of information is available. In particular, depending on 
the sign of the product an, the nonlinearity will either 
counteract {an > 0) or enhance {an < 0) the diffrac- 
tive broadening. Furthermore, the velocity of the particle 
bunch is left unchanged, and no asymmetry is introduced 
on an initially symmetric bunch. Of special interest is 
the case an > 0, when shape-preserving soliton solutions 
are possible as a balance between linear diffraction and 
nonlinear self-focusing effects. 

However, the properties of the full Eq. ^ are not 
known and will be the subject of the present work. In 
order to see the physical significance of the new term, it 
is instructive to qualitatively discuss the nonlinear po- 
tential, V{ip), which in the case a < is given by 




(2) 



Consider first the case — s- and aK > 0, i.e., assume 
that also k < 0. Then, sech-shaped solutions form a 
well-shaped potential that allows bound states, solitons, 
to exist. The fundamental soliton solution corresponding 
to Eq. Ill) is 



= Aq sech(aa;)e 
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The nonlocal part of the potential introduced by /i 7^ 
contributes a monotonous term to the potential and cre- 
ates an asymmetry. A qualitative plot of the total poten- 
tial corresponding to a field shaped as the fundamental 
soliton (choosing Aq = 1, a — I, and n — —1), is shown 
in Fig.n using different values for /i. It is clear that if /i 
is small, the evolution of an initially soliton-shaped pulse 
should involve an acceleration in a direction determined 
by the sign of /i, but the change of the pulse shape can be 
expected to be slow due to the similarities with the con- 
ditions for soliton propagation. For large however, it is 
obvious that there can be no pulse-shaped stationary so- 
lutions, since the total potential then is monotonous, and 
thus is unable to provide the compression effects needed 
for confinement. By noticing that the slope of the po- 
tential varies over the pulse, and that the more intense 
parts are accelerated/decelerated more, strong internal 
pulse dynamics is anticipated. Our subsequent analysis 
will confirm this intuitive picture. 



III. PERTURBED SOLITON DYNAMICS 

From the potential picture it is expected that one of the 
main effects of the nonlocal term is to induce an accelera- 
tion of an initially stationary pulse. A more quantitative 
analysis of this effect can be carried out by investigating 
the adiabatic evolution of the soliton solution, Eq. ||3Jl, 
in the presence of a small but finite value of /i. This can 



FIG. 1: A qualitative plot of the potential, Eq. lj2j, for a 
pulse-shaped field using different values of /i. 



conveniently be done using a direct variational approach, 
see, e.g., @. A suitable trial function is 

VtCx, z) = ^sech[a(x - M)]e*[<^(^"*^)+^1, (4) 

where A(z), O'iz), C'(z), D(z), and M{z) are unknown 
parameter functions to be determined by the variational 
procedure. We emphasise that this ansatz function ne- 
glects any asymmetric pulse deformations and conse- 
quently can model only part of the dynamical evolution. 
Using Ritz optimisation, these parameter functions can 
be determined and the following approximate solution is 
obtained 



Aq sech 



X exp 




where 



2aM^o 2 
X — -z . 



(5) 



(6) 



This solution is consistent with the classical NLS equa- 
tion in the sense that the fundamental soliton, Eq. (j^J, is 
recovered in the limit when — s- 0. The solution in the 
general case describes a soliton being accelerated in the 
original frame of reference. The acceleration, 7, is given 
by 7 = AafiA^/i and the concomitant shift of the group- 
velocity is associated with a frequency shift proportional 
to z. 

In order to check the approximate analytical solutions, 
but also to obtain results in parameter ranges where ana- 
lytical solutions are not available, Eq. ^ has been solved 
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FIG. 2: The numerically obtained dynamics of an initially 
soliton-shaped pulse using a weak nonlocal term. 
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FIG. 3: By increasing the strength of the nonlocal term, 
strong pulse shape distortion is introduced. 



numerically. For this purpose, the standard split-step 
Fourier method for handling the NLSE is modified to in- 
clude the effects of the nonlocal term. An example of 
the dynamics caused by a weak nonlocal term is seen in 
Fig. 121 The initial pulse, which is centred around x = 0, 
is the fundamental soliton, making a comparison with the 
variational result simple. The parameters are a = —I, 
K = —I, fi = 0.1, and Aq = I, and the propagation 
distance is L = 30. The pulse is accelerated towards 
negative x coordinates, but the shape is only weakly dis- 
torted. Thus, for small values of the variational result, 
which predicts the final centre position to be a; = —60, 
describes the propagation in an excellent way. 

However, for larger values of /i, the asymmetric defor- 
mation of the pulse becomes more important and tends to 
leave an extended tail behind the main pulse. A numer- 
ical simulation result of this situation is shown in Fig. |31 
where /i = 1. The propagation distance is L = 10, mak- 



ing the final i-position similar to the example above. It 
is seen that the internal dynamics of the pulse is signif- 
icantly stronger; the peak power decreases, and a wake 
field is developed. Again using the variational result, the 
final position is expected to be x = —66.7, but the de- 
caying peak power makes the potential less steep during 
propagation, which causes the acceleration to decrease. 

Clearly, situations where fi is large can not be analysed 
by a variational approach based on the soliton ansatz. In 
order to study the deformation dynamics in some more 
analytical detail, we will instead use a perturbation anal- 
ysis. 



IV. PERTURBATION ANALYSIS 

Although a general solution of Eq. based on ana- 
lytical methods is not possible, the initial dynamics can 
be described using a perturbation analysis. For this pur- 
pose, Eq. ^ is rewritten as a coupled system in the real 
amplitude. A, and the phase, 9, of according to 



dA^ 

80 
dz 
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dx \ dx 



(7) 
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For an initially unchirped pulse, i.e., 9(x, 0) — 0, the ini- 
tial amplitude modulation first creates a phase modula- 
tion proportional to z, which then generates a subsequent 
change, proportional to z^, of the amplitude modula- 
tion. Let us consider the case of the fundamental soliton, 
Eq. as initial field, since the diffraction and the Kerr 
nonlinearity then balance each other. Thus, as initial 
condition we consider A(x, z — 0) — Aq sech(ax), where a 
is related to Aq according to Eq. JSJ and 0{x, z = 0) = 0. 
Using these in the right hand side of Eq. (jHJ, the ini- 
tial evolution of 9 is obtained, and this solution can then 
be used to find the lowest order modifications of A from 
Eq. CZJ). 

However, as found using variational analysis, the pulse 
evolution is, due to the effects of the nonlocal term, most 
conveniently described in an accelerated coordinate sys- 
tem. The proper value of the acceleration can be taken 
from the previous section, but it is also instructive to 
derive it using an analogy with Ehrenfest's theorem in 
quantum mechanics. It is straightforward to show that 
the motion of the mean position, (x), of the bunch obeys 
the equation of motion 
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dz^ 



where the averaging is defined according to 



(9) 



(10) 



4 



0.8- 



0.6 



0.4- 



0.2 







- - initial 






— numerical 






~ ~ psrturbational 


- 


i / \ \ 






/ / \ \ 
/ / \ \ 





FIG. 4: Comparison of the asymmetric shapes predicted by 
analysis and numerical simulations, respectively. 



and the force is F = —dV/dx. The acceleration obtained 
in this way is identical to that derived using the varia- 
tional approach. Thus, 70 = 7 and the new coordinate, 
^, is defined according to Eq. ©. The amplitude and 
phase are then obtained as 



A = Ao sech(a^) 



X y 1 + -iaa^Al tanh(a^) 



sech^(a^) 



a 



[tanh(aO + 1] ^ ^ 



z2,(ll) 



(12) 



The first part of the phase, Eq. H12|l . does not depend 
on ^, and is identical to the phase of the fundamental 
soliton. The second part is tanh-shaped, which is due 
to the form of the nonlocal potential term. As expected 
from the potential picture, it is found that the amplitude 
becomes asymmetric. 

Due to the limited accuracy, the perturbation anal- 
ysis can only be applied within a certain propagation 
distance, obtained from Eq. I|ll|) as 



AaafiAlz^ < 1 



z <^l/^J\4aanAl\. (13) 



Using the same numerical parameters as above, no signifi- 
cant changes in the amplitude are seen within that range. 
However, as seen in Eq. ()13|l . the application range de- 
creases slowly as fi increases. Thus, a large value, /i = 10, 
has been used in Fig. ^ where the perturbation analysis 
is compared with the numerically obtained result after a 
propagation distance L = 0.2. It is seen that the pulse is 
starting to "lean to the side" , and that the perturbation 
profile is in good agreement with the numerical result, 
although its peak value is slightly too large. 



V. WAVE-BREAKING 

The nonlocal potential term, proportional to fi in 
Eq. (O, gives rise to a force, F(^) ~ — This im- 
plies that the central parts of the bunch are affected 
by a stronger force than the wings, and will acceler- 
ate/decelerate more. In fact, this is the basic mechanism 
behind the steepening and the deformation of the bunch. 
It is clear that after a certain distance of propagation, 
the high amplitude parts should overtake/be overtaken 
by the low amplitude parts of the bunch. However, the 
finite diffraction will prohibit the development of an in- 
finite amplitude gradient, and the "overtaking" between 
different parts of the bunch instead leads the appear- 
ance of oscillations on the amplitude at the base of the 
steepening side of the bunch. This feature is completely 
analogous to the wave-breaking phenomenon in nonlinear 
defocusing Kerr media, 9] , with the difference that in the 
latter case, the corresponding force is an odd function, 
which implies that the wave remains stationary and that 
the wave-breaking phenomenon occurs symmetrically on 
both sides of the pulse. 

In order to estimate the length-scale of the wave- 
breaking phenomenon, the perturbation solution for the 
amplitude, Ec^. Hll|) . can be used. Thus, the order of 
magnitude of the wave-breaking distance, z^b, is esti- 
mated as the shortest propagation length for which the 
amplitude has a zero. This is easily shown to occur at 
Zujf, = •\/3/(4|a|a/i^Q). It is interesting to note that this 
approach gives the same result as the one used in [T?il |. 
which was based on the local velocity shear in the pulse 
created by the nonlinearly induced chirp, provided the 
latter is generalised to include the mean acceleration of 
the pulse. 

By increasing the propagation distance in Fig. 01 oscil- 
lations on the amplitude will start to occur at the base 
of the pulse on the steepening side, i.e., close to x — —3. 
Numerically we define the wave-breaking distance as the 
propagation distance where the amplitude acquires a sec- 
ond maximum. In Fig. [S] the analytical prediction for 
the wave-breaking distance is compared with the result 
of the numerical computations. The results show very 
good agreement, although the numerical results tend to 
be somewhat larger than predicted. However, since only 
an order-of-magnitude estimate has been made, the re- 
sult is quite satisfactory. In particular, the analytic result 
predicts very well how the wave-breaking distance scales 
with fi. 



VI. STATIONARY SOLUTIONS 

As already discussed, the purely reactive case, corre- 
sponding to /i = 0, allows a soliton solution, Eq. ||2Jl, 
containing a finite number of particles. In order to in- 
vestigate whether similar stationary solutions exist also 
in the general case, we return to Eqs. {Tj) and ©, which 
describe the evolution of the (real) amplitude and the 
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FIG. 5: The wave-breaking distance as predicted by analysis 
and numerical simulations, respectively. 
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FIG. 6: A numerically obtained step-like solution. 



phase of the wave, respectively. Based on our previ- 
ous results, we will look for solutions that are stationary 
in an accelerated frame of reference, i.e., we introduce 
^ = x — 7z^/2, where the acceleration, 7, now is unknown 
and has the character of an eigenvalue. Stationarity im- 
plies that the amplitude depends only on the coordinate 
^, i.e., A = A{^). The phase variation can then be found 
explicitly, and the system becomes 



= Co 



a 



(PA 



CiA- 



2a 



12a ' 
-kA^ 



(14) 

HA A'^dC = 0,(15) 



where Ci is a constant, which acts as a second eigenvalue. 
However, by rescaling a, 7, k, and /i, the equation can 
be rewritten as 



fA 



A 



2a 



kA^ + ^lA t A^d£,' = 0. 

J —OC 



(16) 



In effect, this normalisation sets Ci = 1, and the phys- 
ical significance of this can be found by letting /x ^ 
(implying also that 7^0). The fundamental soliton is 
then recovered from Eq. (|15() . and the phase is given by 
Ciz. Thus, setting Ci = 1 corresponds to normalising 
with respect to the propagation constant. 

Assume that a pulse-shaped stationary solution exists. 
This implies that A — > as ^ — > —00, and asymptotically 
the field should satisfy the equation 



PA 

'de 



A 



2a 



0, 



(17) 



which can be solved in terms of the Airy functions, Ai{x) 
and Bi(a;), as 



v4 = Di A 



74 



2a 



\{/2^ 



D2 Bi 



7^ - 2a 



{/2a^^ 



(18) 



In order to obtain a solution containing a finite number of 
particles, it is necessary that 7 < and that D2 — 0. For 
pulse-like solutions, the amplitude of the solution must 
also vanish as — + cxd and the corresponding asymptotic 
equation is 



a^ + il + ,W)A-ll-=0. (19) 

Here, W is the total number of particles, which has been 
assumed to be finite, and the corresponding solution is 



A = D-i Ai 



+D4Bi 



- 2{1 + ^iW)a \ 
'7^-2(l + /^T4^)a' 



(20) 
(21) 



Since the acceleration has already been chosen to be neg- 
ative, this implies that the asymptotic solution will be the 
sum of two oscillating Airy functions as f — > 00. How- 
ever, the total number of particles of such a solution is 
infinite, and a contradiction has been reached. Thus, we 
conclude that there are no stationary solutions to Eq. 
containing a finite number of particles. 

On the other hand, if the condition A as ^ 00 
is relaxed, step-like solutions can be found. By assuming 
in Eq. H16|) that A — > Aoc when ^ — > 00, the integral term 
is asymptotically equal to /iAj^^. This term can cancel 
the term that gives rise to the Airy solutions, provided 
that Aoo — \J^I (2a^). A solution of this type has been 
calculated numerically by using a = — 1,k = — 1,7 = 
— 1, and /i = 1, and by choosing the amplitude for the 
asymptotic solution for negative ^. The result has been 
plotted in Fig. El and it is seen that the predicted value 
for the asymptotic amplitude, A^a = l/v^, is correct. 
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VII. CONCLUSION 

In conclusion, a generalised NLSE, describing the non- 
linear longitudinal dynamics of high energy charged par- 
ticle beams in accelerators within the TWM approach, 
has been analysed using both analytical and numerical 
methods. It has been discussed in qualitative physical 
terms how the inclusion of the resistive part of the cou- 
pling impedance gives rise to both an acceleration and 
a deformation of the particle bunch. These effects have 
been analysed analytically using both a variational anal- 
ysis and a direct perturbation analysis of the initial dy- 



namics, and the results have been shown to be in good 
agreement with numerical simulations. It has also been 
shown that for impedances with a large resistive part, 
the deformation leads to self-steepening and eventually a 
wave-breaking phenomenon similar to that occurring in 
nonlinear optics. The scale length for this effect has been 
estimated and has also been shown to be in good agree- 
ment with numerical results. Finally, it has been shown 
that no stationary pulse-like solutions with a finite num- 
ber of particles exist for the generalised NLS equation, 
but semi-infinite shock solutions are possible. 



[1] R. Fcdole and G. Miele, Nuovo Cimento D 13, 1527 
(1991). 

[2] R. Fedele, G. Miele, L. Palumbo, and V. G. Vaccaro, 

Phys. Lett. A 179, 407 (1993). 
[3] R. Fedele, F. Galluccio, and G. Miele, Phys. Lett. A 185, 

93 (1994). 

[4] D. Anderson, R. Fedele, V. Vaccaro, M. Lisak, 
A. Berntson, and S. Johanson, Phys. Lett. A 258, 244 
(1999). 

[5] R. Fedele and D. Anderson, J. Opt. B: Quantum Semi- 
class. Opt. 2, 207 (2000). 
[6] R. D. Ruth, Proc. of a topical course of the joint US- 



CERN school on particle accelerators, M. Month and S. 

Turner (editors). Springer- Verlag, Berlin (1989). 
[7] J. D. Lawson, The physics of charged-partide beams 

(Clarendon Press, Oxford, 1988). 
[8] D. Anderson, M. Lisak, and A. Berntson, Pramana 57, 

917 (2001). 

[9] W. J. Tomlinson, R. H. Stolen, and C. V. Shank, J. Opt. 
Soc. of Am. B 1, 139 (1984). 
[10] D. Anderson, M. Desaix, M. Lisak, and M. L. Quiroga- 
Teixeiro, J. Opt. Soc. of Am. B 9, 1358 (1992). 



